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The effects of laser-induced prealignment on tlie deflection of paramagnetic molecules by inhomo- 
geneous static magnetic field are studied. Depending on the relevant Hund's coupling case of the 
molecule, two different effects were identified: either suppression of the deflection by laser pulses 
(Hund's coupling case (a) molecules, such as CIO), or a dramatic reconstruction of the broad dis- 
tribution of the scattering angles into several narrow peaks (for Hund's coupling case (b) molecules, 
such as O2 or NH). These findings are important for various applications using molecular guiding, 
focusing and trapping with the help of magnetic fields. 
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I. INTRODUCTION 



Manipulating the translational motion of atoms and 
molecules by means of inhomogeneous external fields has 
been studied intensively for many years. Since the pio- 
neering work of Stern and Gerlach that demonstrated 
quantization of atomic trajectories in inhomogeneous 
magnetic field [1], the dynamics of many other systems 
has been studied both in electric and magnetic fields. 
An important milestone was, for instance, separation of 
molecules in different quantum states in order to create 
a maser, a molecular amplifier of photons [2|- Nowadays, 
the physics of the deflection of atoms and molecules by 
inhomogeneous fields is as hot as ever, including studies 
focused on the motion in the static inhomogeneous elec- 
tric S-Q and magnetic 3, 8, 9] fields, and even laser fields 
[10l - [l3| . In the case of laser dcfiection, some novel appli- 
cations in molecular o ptic s have recently appeared, such 
as molecular lens [ifl, lU and molecular prism [Tj, [T^ . 
The interaction between a molecule and an external field 
depends upon the orientation of the molecule. The 
field-molecule interactions become intensity-dependent 
for strong enough fields due to the field-induced mod- 
ification of the molecular angular motion [IJ, |l5[. It 
was lately shown that the intensity-dependent molecu- 
lar polarizability-anisotropy provides means for tailoring 
the dipole force felt by molecules in the laser field 13]. 
More recently, a method for controlling the scattering of 
molecules in external fields by additional ultrashort laser 
pulses inducing field-free molecular alignment was sug- 
gested fie- 18]. 

In this work we return to the basics, and study the 
prospects of the ultrafast laser control of molecular dc- 
fiection in the Stern-Gerlach (SG) arrangement. It was 
shown in the past that molecular scattering in magnetic 
fields is affected by rotational alignment caused, for ex- 
ample, by collisions in seeded supersonic beams jl9| . 
Here we demonstrate that this process can be efficiently 
and fiexibly controlled by novel ultafast optical tools al- 
lowing for prcshaping the molecular angular distribution 
before the molecules enter the SG apparatus. This can 
be done with the help of numerous recent techniques for 



laser molecular alignment, which use single or multiple 
short laser pulses (transform limited, or shaped) to tem- 
porarily align molecular axes along certain directions (for 
introduction to the rich physics of laser molecular align- 
ment, see, e.g. [Ij, [l5|, [20l - |23| b Short laser pulses ex- 
cite rotational wavepackets, which results in a consider- 
able transient molecular alignment after the laser pulse 
is over, i.e., at field-free conditions. In the present paper, 
we will consider only molecules with a permanent mag- 
netic dipole moment, i.e., open shell molecules. The open 
shell molecules are classified into Hund's coupling cases 
according to their angular momenta coupling [2, |24| . In 
the Hund's coupling case (a), the angular momentum of 
electrons and their spin are coupled to the internuclear 
axis, while in the Hund's coupling case (b), the electronic 
spin and internuclear axis are not strongly coupled. We 
will consider magnetic deflection of different paramag- 
netic molecules subject to a short prealigning laser pulse. 
In the Hund's coupling case (a), the magnetic moment 
is coupled to the internuclear axis, and by rotating the 
molecule, one rotates the magnetic moment as well. This 
reduces substantially the Zeeman effect and effectively 
turns off the interaction between the molecule and mag- 
netic field (like a rotating electric dipole that becomes 
decoupled from a static electric field [131 )• I'^ the Hund's 
coupling case (b), the magnetic moment is barely coupled 
to the internuclear axis. However, laser-induced molecu- 
lar rotation creates an effective magnetic field which adds 
to the SG field and modifies the defiection dynamics. As 
a result, as we show below, a broad and sparse distribu- 
tion of the scattering angles of defiected molecules col- 
lapses into several narrow peaks with controllable posi- 
tions. 

The paper is organized as following. In Sec. |lT] we 
outline the general theorical framefork: first, we briefly 
discuss the SG deflection mechanism (Sec. Ill Ap . and pro- 
vide several needed facts on the laser-induced field-free 
alignment in Sec. IIIBl Then, the interaction details for 
the Hund's coupling case (a) and Hund's coupling case 
(b) are given in Sec. Ill CI and III D[ respectively. Further 
discussion of the Hund's coupling case (b) and hyper- 
fine structure appears in the Appendix in Sec. IVII In 
Sec. mil we apply the above theoretical tools to the laser- 



controlled magnetic scattering of CIO (Sec. IIII Ap . O2 
(Sec. ImBj) and NH (Sec. HITOI) molecules. Discussion 
followed by the conclusions, are presented in Sec. llVl and 
IVl respectively. 



II. GENERAL THEORY 
A. Stern-Gerlach deflection 

Once a beam of molecules enters into a SG magnetic 
field, the initial eigenstates of the system adiabatically 
become \'<i'i{B)): 



where 9 is the angle between the molecular axis and the 
polarization direction of the pulse, a||,a_L are the par- 
allel and perpendicular polarizability components, and 
e is the femtosecond pulse envelope. Since the aligning 
pulse is short compared to the typical periods of molec- 
ular rotation, it may be considered as a delta-pulse. In 
the impulsive approximation, one obtains the following 
relationship between the wavefunction before and after 
the pulse applied at t = (see e.g. [23|, and references 
therein): 



^{t^0+) = exp[iPcos^6']*(t = 0" 
where the kick strength, P is given by: 



(6) 



|*.(i3))=5]a,(i3)|vl',). 



(1) 



where the coefficients aj{B) depend on the magnetic field 
B as a parameter, and |^j) is a basis for the free molecule. 
In this work we consider the magnetic field to be: B — 
B{z)z, i.e., pointing in the z direction, with a practically 
constant gradient along the z direction in the relevant 
interaction region. 

The force acting on the molecule is given by: 

^^^^ (2) 



Force 



-\/E ■ 



dB dz' 

where E is the energy of the molecule. The deriva- 
tive dE/dB may be obtained by means of the Hellman- 
Feynman theorem, that is, for a system being in the i-th 
energy eigenstate of the system, the force is proportional 
to: 



fifj 

= {^^{B)\^\^,{B)) 



(3) 



where H^ is the Zeeman term of the Hamiltonian. Since 
Hz is proportional to B, we conclude that a molecule in 
an energy eigenstate will be deflected by a force that is 
proportional to: 



A^(*.(s)l^l*.(i?)>. 



(4) 



Eq. m will allow us to consider the distribution of 
forces. In order to take into account the absolute amount 
of deflection, though, one has to consider the field gradi- 
ent as well (Eq. [2]). For more details, see, e.g. Q]. 



B. Laser-induced field- free alignment 

If the molecules are subject to a strong linearly polar- 
ized femtosecond laser pulse, the corresponding molecule- 
laser interaction potential is given by: 

Hml = --7e'^[{a\\- a±)cos'^e + a±] , (5) 



P = (1/4/1) • (an - ai_) 



\t)dt. 



(7) 



We assume the vertical polarization of the pulse (along 
the z-axis, and parallel to the SG magnetic field). Phys- 
ically, the dimensionless kick strength P, is equal to the 
typical amount of angular momentum (in the units of K) 
supplied by the pulse to the molecule. In order to find 
'^{t = 0+) for any initial state, we introduce an artificial 
parameter ^ that will be assigned the value ^ = 1 at the 
end of the calculations, and define: 

*5 = exp [(zPcos^ e)(]-^{t = 0-) = ^ c,(0|*.). (8) 

i 

By differentiating both sides of Eq. |8] with respect to 
^, we obtain the following set of differential equations for 
the coefficients &: 



-^^E' 



(*,, 



I*,: 



(9) 



where c — dc/d^. Evaluation of the matrix elements 
in Eq. [9] is easily obtained by means of the relationship: 
cos^ 6 — (2_Dqq-|-1)/3, where D^ is the rotational matrix. 

Since ^^^o = ^{t = 0") and ^^^i = *(i = 0+) (see 
Eq. [SJ, we solve numerically this set of equations from 
^ = to ^ = 1, and find *(t = 0+). It turns out that 
the population of rotational levels of the kicked molecules 
has a maximum at around hP. 

Finally, we derive the distribution of forces acting on 
a thermal ensemble of molecules pre-aligned by a laser 
pulse. For this, we start from a single eigenstate of a free 
system, apply an alignment pulse in the z direction, and 
then adiabatically increase the magnetic field (in order 
to imitate a smooth process of the molecular beam in- 
jection into the SG deflector). The distribution will be 
proportional to: 



/(-4) = E 



exp 






«J 



'Sa.a, 



(10) 



where fcs is the Boltzmann's constant, Qrot is the parti- 
tion function, i denotes the quantum numbers associated 
with the initial eigenstates of free molecules, Cj denotes 
the coefficients of the free eigenstates that were excited 
by the laser pulse applied to the initial eigenstate i, and 
Aj are the associated matrix elements given in Eq. |4] 
(proportional to the force), between the states adiabati- 
cally correlated with the free states j. 



C. Hund's coupling case (a) 

In this subsection we concentrate on the ^^ClO 
molecule, which presents a good example for the Hund's 
coupling case (a). Denoted as ^11 in its electronic ground 
state, it has a nuclear spin / = 3/2, and it was studied 
well in the past [26-28]. 

In the Hund's coupling case (a), the electronic angular 
momentum and spin are strongly coupled to the internu- 
clear axis, and in the case of CZO, its effective Hamilto- 
nian is given by |24| : 



H, 



eff 



— Hr 



H, 



hf 



Hr 



(11) 



where Hrso is the rotation and spin-orbit coupling, Hhf is 
the magnetic hyperfinc interaction, and Hq is the electric 
quadrupole interaction. Here 



the nuclear and electron g factors, respectively, hn and 
fiB are the nuclear and electron Bohr magnetons, respec- 
tively, and /Uo is the vacuum permeability. 

All the matrix elements for the Hund's coupling case 
(a), including those for the quadrupole interaction, are 
given in |23|. The CIO constants were taken from \24 . 

PT^ I L 

mm- 

When considering the Zeeman Hamiltonian, we will 
concentrate only on the two major terms related with 
electronic angular momentum and spin: 

Hz = fiBT\B) ■ T\L) + gsf^BT\B) ■ T\S). (15) 

The corresponding matrix elements (see Eq. |4]) are given 
in [24]. 

In order to consider the effect of laser-induced align- 
ment (see Sec. IIIBl Eq. [3]), we have derived the following 
matrix elements: 



( v,A;S,J:;J,n,I,F,Mp\Dl*\ri,A;S,J:;J',n,I,F',Mp) 

_ / -.nF-A/f f P 2 F' \ . ,^F' + ,/+/+2 

- ^ ^) y -Mf Alp J^ ' 

X ^{2F' + l){2F+l)i-^p ^j 2|(-1)''-^' 



J 2 J' 
-f) 17 



") V(2J+1)(2J' + 1). 



(16) 



Hrso - Br {T\J) - T\L) - T\S)Y+AT\L) ■ T\S), 

(12) 
where T-^{) is a spherical tensor of rank 1, Br is the ro- 
tational constant in the lowest vibrational level, and A is 
the spin-orbit coupling constant. L and S are the elec- 
tronic angular momentum and spin operators, respec- 
tively. The total angular momentum isJ = N-|-L-|-S, 
where N is the nuclei angular momentum operator. The 
Hund's coupling case (a) basis looks like this: 



\r,,A;S,i:;J,n,I,F,MF) 



(13) 



where rj represents some additional electronic and vibra- 
tional quantum numbers, S and A are the projections of 
the electronic spin and angular momentum on the inter- 
nuclear axis, respectively. For CIO molecule, S = 1/2, so 
that E = ±1/2 and A = 1. The quantity O is 17 ee E + A 
(n = 3/2,1/2, the 3/2-state has a lower energy), and 
F = J + I. 

The Hhf Hamiltonian is given by: 

Hhf = HiL + Hp + Hdip 

= aT^(I)-r^(L) + 6^^1(1) •Ti(S) 

- yiOgsA'B3JvMiv(/^o/4^)ri(I)-ri(S,C2). (14) 

The first term represents the orbital interaction, the 
second one accounts for the Fermi contact interaction, 
and the third term describes the dipolar hyperfine in- 
teraction. Here a and hp are constants, g^ and gs are 



D. Hund's coupling case (b) 

We will continue by discussing the Hund's coupling 
case (b), and consider the Oxygen molecule, in its pre- 
dominant isotopomer ^^O^^O. This molecule is proba- 
bly the most important species among '^E ground state 
molecules, and it was one of the first molecules studied 
in detail [a, l29|. It is a homonuclear diatomic molecule, 
where only odd N's appear because of the Pauli's princi- 
ple and symmetry. This molecule is described well by the 
Hund's coupling case (b), with the effective Hamiltonian 



^eff — ^rot + Hss + Hsr- (17) 

Let us describe separately each term in Eq. [T71 Here 



Hrot = Br'N^ - DN'^ 



(18) 



is the energy of the nuclei rotation, where D is the cen- 
trifugal distortion coefficient. In addition, 



Hss = -5.VI(Mo/47r)V6r2(C) • T^S,, S2) 



(19) 



is the electornic spin-spin dipolar interaction. T^() is a 
spherical tensor of rank 2. 81,82 are electronic spin oper- 
ators. T^{C) = {C^{e,<l))R-^), where C^ is the spherical 
harmonics, and 



Hsr^7T\N)-THS) 



(20) 



is the electronic-spin rotation interaction. The Hund's 
coupling case (b) basis looks like this: 



\r,,A;N,A;N,S,J,Mj). 



(21) 



Here N is the nuclei rotational quantum number, A — 
in our case, S is the electronic spin, which is 1 in our 
case, J — N+S, and Mj is the projection of J onto a fixed 
z-direction in space. All the needed matrix elements and 
constants are given in 24|. The Zeeman Hamiltonian is 
given by: 



Its matrix elements (Eq. H]) are given by: 



(22) 



d{ ry. A; N, A; N, S, J, Mj\T^{S)\r^, A; N' , A; N' , S, J', Mj) 



d(-l) 



J-Mj 



J 1 J' 

-Mj Mj 



JN,N 



(-1) 



J'+S+l + N 



X V(2J' + 1)(2J+1) 
X ^S{S + l)i2S+l), 



S J' N 
J S 1 



(23) 



where d = gsfJ-B- Finally, in order to account for the 
laser- induced prealignment, we derived the following re- 
lation (to be used in Eq. [9]): 



{r^,A;N,A;N,S,J,Mj\Dl*\v,^;N',A;N',S',J',Mj) 



(-1) 



J-Mj 



X V(2^' + 1)(2J+1) 

A^ 2 iV' 



J 2 J' 

Mj Mj 

N' J' S 
J N 2 



(-i)^-A 



^A A 



v/(2iV+l)(27V' + l). 



(24) 



In the case of the oxygen molecule, there is a rela- 
tively strong effect of the spin-spin interaction, which 
complicates our analysis. Therefore, we have also cho- 
sen an additional "^S molecule, ^'^NH for our study. For 
this molecule the ratio between the spin-spin and spin- 
rotation interactions is reduced (compared to the O2 
case) . This makes NH a simpler candidate to test our ro- 
tational effects. The NH molecule was thoroughly stud- 
ied in the past p30l - l33| . and its effective Hamiltonian is: 



H, 



eff 



Hr 



Hss + Hsr + HhFS-> 



(25) 



where Hrot, Hss and Hsr were defined in Eq. [T8lfT9l and 
[ini Since NH has non-zero nuclei spin {N has nuclear 
spin I = I, H has / — 1/2), then it has a hyperfine 
structure described by the Hamiltonian Hhfs- Further 
elaboration on the hyperfine structure of NH (including 
details on the Zeeman term, and the matrix elements 
related to Eq. [S]), is given in the Appendix in Sec. I VII 



III. LASER CONTROL OF THE 
STERN-GERLACH SCATTERING 

A. CIO 

In this part of the work, we apply the theoretical tools 
that were presented in the previous sections to the SG 
scattering of the CIO molecule. This molecule exhibits 
a good Hund's coupling case (a), and details about it 
were already given in Sec. Ill CI We will consider here its 
ground state (T = Oi^), for which A = l,I] = l/2,fi = 
3/2, J = 3/2, F ^0, Mf = 0. In Fig. [U we present the 
force distribution (Eq. [TU| for a CIO molecule in the 
ground state that is deflected by a SG magnetic field. As 
only single molecular state is occupied, the force has a 
well defined single value. 
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FIG. 1: The force distribution for a beam of CIO molecules 
that are deflected by a magnetic field of IT. The temperature 
is OK, therefore only the ground state is considered, and the 
distribution reduces to a single- value peak. 

As the next step, we assume that the molecules are 
subject to a short laser pulse with a kick strength of 
P = 30 (Eq. [7]) before they enter the SG magnetic field. 
The new force distribution is given in Fig. [51 
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FIG. 2: The force distribution for a beam of prealigned CIO 
molecules. The temperature is OK, and the kick strength of 
the laser is P = 30. The prealigned molecule is deflected by 
a magnetic field of IT. This distribution should be compared 
to the one from Fig. [T] 

By comparing Fig. [T]to Fig. [31 it can be observed that 
the effect of the laser-induced field-free alignment is to 



effectively turn-off tlie interaction between tlie molecule 
and the magnetic field. This effect is similar to the one 
discussed by us recently in connection with the scattering 
of polar molecules by inhomogeneous static electric fields 
[17| . Moreover, rotation- induced dispersion in molecular 
scattering by static electric fields was used as a selection 
tool in recent experiments on laser molecular alignment 
[3] . A related phenomenon of the reduction of the electric 
dipole interaction in highly excited stationary molecular 
rotational states was observed there. Further details and 
discussion about the CIO magnetic deflection is provided 
in Sec. HVl 
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B. O2 

In this sub-section we consider the O2 molecule. The 
O2 molecule is described well by Hund's coupling case (b) 
scheme, and the details about it were given in Sec. IIIDI 
First, we consider a beam of O2 molecules at OK, i.e., in 
the ground state {N — 1, J — 0, and Mj — 0). These 
molecules enter a magnetic field of IT, and are deflected 
by this field. The force distribution for these molecules 
is given in Fig. |31 
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FIG. 4: The force distribution for a beam of prealigned O2 
molecules. Different kick strengths (P = 10, 30, 70) are con- 
sidered and the magnetic field is IT (temperature is OK). As 
the excitation is increased, the two major peaks become closer 
to each other. This distribution should be compared to the 
one from Fig. [S] 
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FIG. 5: The force distribution for O2 molecules. The temper- 
ature is 5K and the magnetic field is IT. 



FIG. 3: The force distribution for a beam of O2 molecules 
that are deflected by a IT magnetic field. The temperature 
is OK, i.e., only the ground state is populated and therefore 
the distribution reduces to a single-value peak. 

Second, we consider the action of the prealignment 
pulses of different kick strengths (P = 10, 30, 70) before 
the molecules enter the deflecting field. The distribution 
of forces at OK is given in Fig. 21 where two major peaks 
are observed. As the strength of the pulses is increased, 
higher rotational states are excited, and the peaks be- 
come closer to each other. 

Third, we consider deflection of thermal molecules 
without and with prealignment, in Figs [5] and [HI respec- 
tively. By comparing Fig. |3] to Fig. [BJ we observe an 
additional peak in Fig. [S) As the strength of the pre- 
alignment pulses is increased, the peaks in Fig. |6] are 
changed: they become narrower and the two left peaks 
become closer to each other. Further discussion on O2 
will be provided in Sec. IIVI 



NH 



Finally, we consider the NH molecule. This molecule 
is described well by a Hund's coupling case (b) scheme, 
similar to the O2 molecule, however it has a reduced 
spin-spin to spin-rotation interaction ratio. This makes 
the NH molecule a simpler candidate for the theoretical 
analysis. In Fig. [7] we plot the force distribution for the 
ground state iV = 0, J = 1, Fi = 3/2, F = 1/2 molecules 
that were prealigned by laser pulses of different intensity. 
Mp was taken to be 1/2 for certainty, and the case of 
Mp = —1/2 may be considered similarly (with similar 
consequences, as will be described in Sec. ITV)) . 

One may observe the presence of three major peaks 
now (for the O2 molecules in the ground state there were 
only two peaks). As the strength of the prealignment 
pulse is increased, the major peaks are shifted in position. 
An additional difference between the NH and the O2 
molecules is that now the peak to the right is also shifted 
with increasing the strength of the prealignment laser 
pulse. Further elaboration about this molecule, and the 
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FIG. 6: The force distribution for a beam of O2 molecules, 
prealigned by a laser field (P — 10, 30, 70). The temperature 
is 5K, and the magnetic field is IT. Here we observe three 
major peaks. As the laser excitation strength is increased, 
the peaks become narrower, and the two left peaks become 
closer to each other. 
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FIG. 7: The distribution of forces for a beam of NH 
molecules, that were prealigned (starting from the lowest state 
iV = 0, J = 1, Fi = 3/2, F = 1/2) by means of laser pulses 
of different strengths: P = 10 (green), P = 30 (blue) and 
P = 40 (red). Only Mf = 1/2 is considered, at a 2T mag- 
netic field. 



difference between it and O2, is given in Sec. |IV] 



O2 



In Fig. |8]we plot the forces vs. magnetic field, for sev- 
eral values of J. First, we observe that for a high mag- 
netic field all the curves are separated to form a three 
SG splitting pattern Q. In the limit of the low mag- 
netic field (and slow rotations), the energy spectrum of 
the molecule is rather complex due to the spin-spin inter- 
action [2^. At around IT, though, we are in the regime 
where the spin-rotation (Hsr) interaction has a rather 
strong dynamic control: as N is increased (by the means 
of prealignment, for instance) then a sizable shift of the 
force magnitude is observed. This is the origin of the 
behavior of the distributions of Fig. U and Fig. [B] It 
can also be observed in Fig. [8] that the spin-spin term is 
larger than the spin-rotation one, and it shifts the curve 
ioi J = N from the two other curves. We find also that in 
the case of J = A^ + 1,A^— 1, the forces are more suscepti- 
ble to different A^s, which is observed in the distribution 
of forces in Fig. |3]and Fig. [51 
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FIG. 8: Forces vs. magnetic field for the O2 molecule. The 
y axis is given in arbitrary units, the x axis is given in the 
units of Tesla. Blue, green, and red (solid lines) correspond 
to A = 31, J = 30,31,32 {Mj = 0), respectively. Blue, 
green, and red (dashed lines) correspond to A = 71 and J — 
70, 71, 72 (Mj = 0), respectively. The effects of the spin-spin 
interaction reveal themselves in the fact that the upper level 
J = A is well separated from two almost degenerate levels 
with J = A + 1,A — 1. Magnetic field near IT is optimal 
for observing the sensitivity of the defiecting force to the A 
variation. 



IV. DISCUSSION 



CIO 



First we will discuss the CIO molecule, which exhibits 
a good Hund's coupling case (a). Having both electronic 
angular momentum and spin coupled to the intcrnuclear 
axis, rotation of the molecule by means of short laser 
pulses leads to the rotation of the molecular magnetic 
moment as well. The interaction between the SG mag- 
netic field and the rapidly rotating magnetic moment of 
the molecule will be thus averaged to zero, leading to the 
negligible magnetic forces. 



Fig. [5] also allows us to understand the position of the 
peaks in Fig. 2] and Fig. [S] The right peak that appears 
in Fig. inland does not appear in Fig. ^corresponds to the 
J = N states. Due to selection rules (Eq. [24]), the odd 
Js, i.e., the J = N states, are never excited (if we start 
from J = 0, and Mj = at OK). This is why we observe 
only two peaks, i.e., the J = A^± 1 peaks, in Fig. JH Con- 
sidering a deflection of the molecules in the ground state 
alone is important experimentally. Even if one considers 
an experiment atT = 1K (ksT = 20837MHz), then the 
difference between (A^ = 1, J = 0) and (A^ = 1, J = 2) 
(the next energy) is 624:86MHz, which is large enough. 
Though, for IK we should expect a small peak in the 
distribution of forces for J = N states. In the case of 
larger temperature (Fig. |B]), we start from different Ms, 



and also the odd Js are present, therefore, we observe 
the right peak at Fig [5] As the preahgnment becomes 
stronger, the distribution transforms into three peaks, 
each correspond to either J = N, J = N—1 or J = N + 1 
states. 



C. NH (and the imaginary O2 molecule!) 

Before we start with the NH molecule, we consider an 
imaginary O2 molecule (!). This molecule is similar to 
the O2 molecule, only with a spin-spin interaction that 
is reduced by a factor of 100. The forces vs. magnetic 
field for the imaginary O2 molecule are plotted in Fig. [HI 
where we get a symmetric splitting of the N level into 
J = N,N ±1 {J = N is in the middle, as is intuitively 
expected for SG splittings, unlike in the O2 case). Such 
behavior corresponds to a spin-spin interaction that is 
negligible compared with the spin-rotation. 



3 

O 




0.2 0.4 0.6 O.i 

Magnetic Field [T] 



FIG. 9: Forces vs. magnetic field for the imaginary O2 
molecule (details in the text). The y axis is given in arbitrary 
units, the x axis is given in the units of Tesla. Blue, green, 
and red (solid lines) correspond to A'^ = 31, J — 30,31,32 
{Mj = 0), respectively. Blue, green, and red (dashed lines) 
correspond to iV = 71 and J = 70,71,72 {Mj = 0), respec- 
tively. At about IT for this molecule, we observe approxi- 
mately a symmetric splitting into three graphs for J — N 
and A'^ ± 1, where J = A'^ is in the middle (unlike in the O2 
case) . 

As we intuitively suggested in Sec. HI when one ap- 
plies preahgnment to molecules belonging to the Hund's 
coupling case (b), the electronic spin feels the SG field 
combined with the effective magnetic field due to nuclei 
rotations. The latter field is along the A^-vector, i.e. per- 
pendicular to the molecular axis. A strong enough verti- 
cally polarized laser pulse excites molecular rotations in 
the vertical planes containing the z-axis. As a result, the 
rotation-induced effective magnetic field is perpendicular 
to the vertical SG field. Therefore, the force felt by the 
molecules is given by 



where Kq and Ki are constants (the latter is proportional 
to N or P). 

Fig. |9] presents results of the exact quantum- 
mechanical calculation of the SG force for our imaginary 
O2 molecule. We considered the upper curves in this 
figure, and tried to fit them to the above analytical ex- 
pression. We find an excellent agreement between the 
original data and the fitted curves, and the results of the 
fit are Ki = 0.61,0.25 for N = 71,31, respectively We 
also find a good agreement between the ratio oi N's (i.e., 
71/31 = 2.3) and oi K[s (i.e., 0.61/0.25 = 2.4). 

As we have mentioned before, the NH molecule is also 
characterized by a reduced value of the spin-spin inter- 
action compared to the spin-rotation interaction. There- 
fore, its dynamics should be closer to the imaginary O2 
molecule than to the real O2 molecule considered above. 
In Fig. [TU] we plot some forces vs. the magnetic field 
for the NH molecule, and indeed, observe a triplet-like 
structure similar to that one in Fig. [3] (with the curve 
for J = TV being in the middle for large enough N). 
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FIG. 10: Forces vs magnetic field for the NH molecule (in- 
cluding the fine and hyperfine details), for Ai' = 10 (green), 
N = 30 (blue) and A^ = 41 (red). Only Mf = 1/2 is 
considered here, but Mp = —1/2 gives the same results 
(only higher Alp's will modify the spectrum). For A'' = 10 
the upper/middle/lower curves correspond to J — N,J = 
N — 1, J = A^-l-l, respectively, as in the Oxygen case. For 
N — 30, 41 the upper/middle/lower curves correspond to 
J^N-1,J = N,J^N + 1, respectively, as in Fig. [O] 



By analyzing the results shown in Fig. [TU]we conclude 
that the hyperfine structure details in NH are not cru- 
cially important for our considerations, but the reduced 
spin-spin to spin-rotation interaction ratio for NH de- 
fines the major difference of the deflection dynamics com- 
pared to the case of O2 ■ Also one notices the scaling with 
the magnetic field: here the higher values of the magnetic 
field are required {2T) to observe the collapse of the broad 
distribution of forces into three narrow groups. This is 
due to the increased spin-rotation interaction for NH (as 
compared to 02)- Fig. [TOl explains the behavior of the 
distribution in Fig. [71 where we have noticed that the 
three peaks are shifted as one increases the laser pulse 
strength. 



V. CONCLUSIONS 

We considered scattering of paramagnetic molecules 
by inhomogeneous magnetic field in a Stern-Gerlach- 
type experiment. We showed that by prealigning the 
molecules before they interact with the magnetic field, 
one obtains an efficient control over the scattering pro- 
cess. Two qualitatively different effects were found, de- 
pending on the Hund's coupling case of the molecule. For 
molecules that belong to the Hund's coupling case (a), 
we showed that the deflection process may be strongly 
suppressed by laser pulses. This may be implemented 
as an optical switch in the molecular magnetic decelera- 
tion techniques ^S]. Furthermore, for the Hund's coupling 
case (b) molecules, a sparse distribution of the scatter- 
ing angles is transformed into a distribution with sev- 
eral compact deflection peaks having controllable posi- 
tions. Each peak corresponds to a scattered molecular 
sub-beam with increased brightness. The molecular de- 
flection is considered as a promising route to the sep- 
aration of molecular mixtures. Narrowing and displac- 
ing scattering peaks may substantially increase the effi- 
ciency of separating multi-component beams, especially 
when the prealignment is applied selectively to certain 
molecular species, such as specific isotopes [3J], or nu- 
clear spin isomers [35l . |36| . One may envision more so- 
phisticated schemes for controlling molecular scattering, 
which involve multiple pulses with variable polarization 
for preshaping molecular angular distribution. In par- 
ticular, molecular rotation may be confined to a certain 
plane by using the optical molecular centrifuge approach 
[371 i38| , double-pulse ignited molecular propeller ^ff] , or 
permanent planar alignment induced by a pair of delayed 
perpendicularly polarized short laser pulses [40, |41| . If 
molecules are prepared like this, a narrow angular peak is 
expected in their scattering distribution from a magnetic 
field. The position of the peak is controllable by inclina- 
tion of the plane of rotation with respect to the deflect- 
ing fleld, similar to a related effect for molecular scatter- 
ing in inhomogeneous electric fields see |18|). Moreover, 
further manipulations of the deflection process may be 
considered, e.g., by using several SG fields with varying 
directions. Magnetic deflection of O2 molecules subject 
to laser-induced fleld-free manipulations, is currently a 
subject of an ongoing experimental effort. 
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VI. APPENDIX: NH (HUND'S COUPLING 
CASE (B)) HYPERFINE STRUCTURE 

The hyperflne structure for the NH molecule is de- 
scribed by the following Hamiltonian [2J, ISlMSSj : 



k 

- J2tkVT0T\lk)-T'{S,C'{uj)) 

k 

- eT\S/E2) ■ T\Q^) 

+ J2cjik)T\N)-T\l,), 



(27) 



where the sum over k = 1,2 represents the terms for both 
nuclei. The flrst term is the Fermi contact interaction, 
the second term is the dipolar interaction, the third one is 
the quadrupole term (this term exists only for the ^"^iV), 
and the last term accounts for the nuclei spin-rotation 
interaction. In the calculation of matrix elements we first 
coupled J = S -I- N, Fi — In + 3 and only then F = 
Iat -I- Fi. All the matrix elements are diagonal in F, and 
the first three terms are given in [24j . 

The nuclear spin-rotation interactions are given by: 

(7?, A, N, S, J, h,F^,h,F, Mf\T\'H) ■ T\li) 
\i^,K,N',S,J',h,F[,h,F,Mp) 



(-1) 



r+Fi+h 



h J' Fi 



X VA(/i + I)(2/i + l)<5Ar,Ar,(-I)'^'+^+l+^ 

X ^{2J+l){2J' + l)l^j ^ \ 



X yjN{N+\){2N + l), 
where /i = Ih, and 

{r,,A,N,S,J,h,Ful2,F,MF\T\N)-T\l2) 
\7j,A,N',S,.r,h,Fi,l2,F,MF) 

- (-1)^1+^+^^ / ^2 F{ F]^ 
^ ^ ' [Fi I2 I j 

(2Fi + I)(2F^ + l)|^j ^j ^1^ 

' N' J' S 
J N 1 



(28) 



X (-1)^ +^+^+^^{2J + 1)(2J' + I) 
X SN,N'VNiN + l){2N + l), 



(29) 



where I2 = In- The constants were taken from [32|. In 
the Zeeman Hamiltonian we consider only the contribu- 
tion due to electronic spin: 



Hz^fiB9sT\B)-T\S), 



(30) 



and we neglect other small contributions coming from 
the nuclei's rotation and spin, and electronic anisotropy. 
The Zeeman matrix element is proportional to: 

(7?, A, N, S, J, /i, Fi, /2, F, Mf\T^{S) 
|r,,A,iV,5,J',/i,Fi',/2,F',Mf) 
F 1 F' 
-Mf Mp 






X (-!)"■ +"^+^+^^ ^{2F + 1)(2F' + 1) I 5 ^1 t 

5 J' ^ 
J S* 1 



X (-l)'^'+^+i+^v/(2^' + l)(2^+l) ' ^ "^' ^ 



X \/S(5 + l)(2^+l), (31) 

where it is no more diagonal in F. Finally, for the align- 



ment calculations the following matrix element is useful: 



r;,A,7V',5,J',/i,Fi',/2,F',Mf) 

F 2 F' 
-Mp Mp 



-1 



,F-Mf 



-1 



,F' + Fi+2+/2 



v/(2F+l)(2F' + l) 



F{ F' /2 
F Fi 2 



-l)'^'+^+^+V(2J+l)(2J' + l)H' }J f 
N f N 2 N'' 



)^ ( ^ T ) V(2A^ + l)(2iV'+l)- (32) 
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